Transverse modulational instability of partially incoherent soliton stripes 
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Based on the Wigner distribution approach, an analysis of the effect of partial incoherence on 
the transverse instability of soliton structures in nonlinear Kerr media is presented, ft is explicitly 
shown, that for a Lorentzian incoherence spectrum the partial incoherence gives rise to a damping 
which counteracts, and tends to suppress, the transverse instability growth. However, the general 
picture is more complicated and it is shown that the effect of the partial incoherence depends 
crucially on the form of the incoherence spectrum. In fact, for spectra with finite rms-width, the 
partial incoherence may even increase both the growth rate and the range of unstable, transverse 
wave numbers. 

PACS numbers: 42.25.Kb, 42.65.Jx, 



I. INTRODUCTION 

Nonlinear phenomena like self-focusing, collapse, mod- 
ulational and transverse instabilities of cylindrical light 
beams are some of the most fundamental consequences 
of the interplay between linear diffraction and self-phase 
modulation in nonlinear Kerr media. Various physical 
mechanisms, which tend to suppress such instabilities 
e. g. nonlinear saturation, have been diligently analyzed 
in a number of works, see references in These fun- 
damental instability problems have continued to attract 
attention in connection with new scientific and techni- 
cal developments. There is currently a strong interest 
focused on the effects of partial incoherence on different 
nonlinear instabilities 0. y ■ 0- [!■ The results of these 
studies show that the modulational and collapse instabil- 
ities tend to be suppressed when the waves are partially 
incoherent. Recently, the effect of partial incoherence on 
the transverse modulational instability of soliton stripes 
in nonlinear Kerr media has been investigated, see 
A soliton stripe is a semi-localized structure, which is of 
self-trapped soliton form in the x direction, uniform in 
the y direction and propagates in the z direction. While a 
ID soliton is resilient to perturbations, the soliton stripe 
exhibits instability with respect to transverse perturba- 
tions, i. e. perturbations in the y direction, see e. g. 0. It 
has been shown, that when the stripe is partially in- 

coherent in the y direction, the transverse modulational 
instability tends to be suppressed and the break-up of 
the stripe, due to the transverse modulational instabil- 
ity, can be prevented provided the incoherence is suffi- 
ciently strong. This behavior is similar to that of the ID 
modulational instability. However, analysis of the trans- 
verse modulational instability is more complicated than 
the corresponding analysis in the case of ID modulational 
instability. In fact, even in the fully coherent problem, 
the problem of finding the growth rate as a function of 
the wave number of the perturbations does not have an 



explicit analytical solution, cf |3- 

In the present work we present an analytical investiga- 
tion of the effect of partial incoherence on the transverse 
instability of soliton structures in nonlinear Kerr media. 
It will be shown that in the case of a Lorentzian incoher- 
ence profile, the growth rate of the transverse instability 
can be expressed simply as the growth rate for the co- 
herent case minus a stabilizing damping rate due to the 
partial incoherence. However, we also show that the case 
of a Lorentzian profile represents a very special case and 
the effect on the growth rate in a general case depends 
crucially on the form of the incoherence spectrum 0. 
Using a perturbation approach to the dispersion relation 
for a general form of the incoherence spectrum, we show 
analytically that for weak incoherence spectra of finite 
rms-width, the region of instability always widens and 
the growth rate is increased in some part of the region. 
This result agrees well with a recent numerical study of 
the transverse instability of partially incoherent solitons 
H , where the angular spectrum is assumed to have Gaus- 
sian form. 



II. THE WIGNER APPROACH 

Our analysis is based on the Wigner approach, which 
has been shown to be a convenient tool for analyz- 
ingthe dynamics of partially incoherent light waves, cf 
.1, 4, 5. 6j. 

The starting point of the analysis is the nonlinear 
Schrodinger (NLS) equation for the complex wave field, 
i/j(r,z), describing the two dimensional propagation of 
a partially coherent wave in a diffractive nonlinear Kerr 
medium, 
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where the bracket (•) denotes statistical average, z is the 
distance of propagation and r = (x, y) denotes the trans- 
verse coordinates. The medium response is here assumed 
to depend only on the statistically averaged intensity i. e. 
I = (ijjip*). This form of the NLS equation is valid when 
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the medium response time is much larger than the char- 
acteristic time of the stochastic intensity fluctuations and 
yet much shorter than the characteristic time of the wave 
envelope variation. 

Within the Wigner approach, Eq.JTJ is transformed into 
the Wigner-Moyal equation for the corresponding Wigner 
function p(r,p, z), viz 

|+P"|+2iV(r, 2 )s in ^.|)p(r,p ) ,)=0. (2) 

The Wigner distribution is determined by the stochastic 
properties of ip(r,z) and conversely N = (\ip\ 2 ) = J pdp 
is the average field intensity. In the present application 
we consider a background solution in the form of a soliton 
stripe, i. e. a semi-localized structure, which constitutes 
a self-trapped soliton form in the x direction, is uniform 
in the y direction and propagates in the z direction. This 
structure is assumed partially incoherent in the y direc- 
tion. The corresponding intensity and the concomitant 
Wigner distribution are 



N (x) = sech 2 (x) 



(3) 



and 



Po(x,p) 



2 sm(2xp x ) 
sinh(2x) smh(irp x 



-G(p y ) = R (x,p x )G(p y ), 



(4) 

respectively, where G{p v ) characterizes the spectrum of 
the partial incoherence in the transverse direction. In or- 
der to analyze the stability of this background solution, 
we consider the dynamics of a small perturbation by writ- 
ing p = po(x, p) + pi(r, p, z), where p\ <C p . The linear 
evolution of the small perturbation p\ is then governed 
by the equation 

dpi dpi _ . ( I d d \ 

^ +P . — + 2iV sm^-— jp l( r,p,,) + 

/ 1 ~d\ 

2ni(r,z)sin I -— ■ — I p (x,p) = 0. (5) 

where n\ = J pidp. When considering the transverse 
modulational instability, the perturbations can be as- 
sumed to be described by harmonic variations, i. e. 
ni(r, z) — n(x) cos(fcy) exp(Tz), where k is the wave 
number of the transverse perturbation. With this ansatz 
for the perturbation, Eq.Q can be rewritten in the com- 
pact form 

Q PJL + p . ^1 + 2No s pi + G+ tiSRq + ^^CR = 
oz or k oy 

(6) 

where we have introduced the operators 



S 



1 d d 

2 dx dp x 



C 



1 d d 

2 dx dp x 



and used the notations G+ = G(p y + k/2) + G(p y — k/2) 

and G- = G(p y + k/2) - G(p y - k/2). 

The solution of Eq. © can be represented as 

pi = [J7(a;,p)cos(fcy) + V r (a;,p)sin(fcj/)]exp(rz), (7) 

where the unknown functions U and V satisfy the equa- 
tions 



TU + p x ^- + 2N SU + kpyV = -G+ nSR 0l 
Ox 

dV 

TV+p x — + 2N SV - kp y U = G- nCR a . 
Ox 



(8) 



This equation system has to be solved subject to 
the consistency conditions JUdp x dp y — n(x) and 
/ V dp x dp y = 0. 



III. THE CASE OF LORENTZIAN 
INCOHERENCE SPECTRUM 

For the development of our analysis it is useful to first 
reconsider the case of a fully coherent wave. The trans- 
verse coherence spectrum is then a Dirac delta function 
i. e. G(p y ) = S(p y ). The earlier introduced notations G+ 
and G_ now become a sum and a difference, respectively, 
of two translated delta functions. The p y dependence of 
the U and V functions can be expressed in similar man- 
ner, i. e. U = [S(p y + k/2) + 5(p y — k/2)]u(r,p x ) and 
V = [5{p y + k/2) - 5{p y - k/2)]v(r,p x ). The combina- 
tions of delta functions now appearing in Eq.JSJ) can be 
shown to be separable, and the resulting system of equa- 
tions reduces to 



Lu v = —nSRn 

2 



Lv + —u = nCRo 



(9) 



where we have introduced a p^-independent operator L 
defined as L = T + p x d/dx + 2N S. Eqs.Q can be 
combined into a single equation for u, which reads 

u - ^p-^LjnCRa + 2Lp- 1 {L}nSR = 0. (10) 

where P _1 denotes the inverse of the operator P{L} = 
k 4 /4 + (L) 2 and curly brackets denote the argument 
of the operator. The solution of the eigenvalue prob- 
lem (r = r(0, k)) cannot be found analytically, and re- 
sort must be taken to approximate analytical techniques 
and/or numerical computations, cf |9l Il0|. As an exam- 
ple, a derivation inspired by direct variational methods 
is given in the Appendix. 

With this result in mind for later comparison, we turn 
back to the partially incoherent problem. In the same 
way as for the coherent case, we can eliminate the func- 
tion V in Eq.© to obtain 

L 2 U + kpy{kp y U + G- nCR ) = -LG+ nSR . (11) 
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Integrating this equation over p y -space we obtain: 



U dp y 



A G_ kpyuCRo dp y 



(12) 



A^G+LnSRo^py = 0, 



where yet another new operator, A, has been intro- 
duced. It is defined as A — [k 2 p 2 + L 2 ] and A^ 1 de- 
notes its inverse. Since the p y dependence in the opera- 
tor A^ 1 is multiplicative, some important simplifications 
can be made. For instance, in the second integral of 
Ea. (|12[) . the ordering of the terms may be interchanged 
as A~ 1 G- py-nCRo = G- p y A~ 1 nCRo- For the subse- 
quent analysis we need the eigenvalue rather than the 
operator itself since A~ 1 nGRo = J2 m a m 1 (Py) c rnnCRo. 
The eigenvalue a" 1 corresponding to the operator A" 1 is 
given by a" 1 = l/(A 2 n + k 2 p 2 ) where, in turn, A^ is the 
eigenvalue of the L 2 operator. 

We will now assume the incoherence spectrum to have 
a Lorentzian profile, G(p y ) = po/[^(j>l + Po)]> with the 
characteristic width pg. This assumption has the impor- 
tant consequence that the integrals appearing in Ea. i|12|) 
can be evaluated explicitly to yield: 



a m lG -Py d Py 



/OO 
a^G+dpy 
-OO 



-1 



(fc/2) 2 + (p + \ m /kf 
2(A m + kp ) 



(13) 



k 2 X r 



(fc/2) 2 + (p + X rn /kf 



(14) 



Thus the dispersion equation, Ea. H12fl . can be expressed 
in the following form 



U dp y -k 2 p- l {kp + L}nCR 

+ 2 (kp + L] P^ikpo + L}nSR = 0. 



(15) 

A comparison of the two dispersion relations (the one 
for the coherent, Ea. (|10f> . and the one for the partially 
incoherent, Ea. (|15fl . case) shows that the only differ- 
ence between the two is the shift in the argument of the 
P operator; the argument L is replaced by (L + kpo) 
in the partially incoherent case. Equivalently, since 
L = T + p x d/dx + 2N S, this implies that T(0, k) = 
r(po,fc) + kpo, where T(po,k) denotes the growth rate 
of the partially incoherent case. Thus we finally come 
to the important conclusion that, for Lorentzian shaped 
incoherence spectrum, the role of the partial incoherence 
on the transverse modulational instability of a soliton 
stripe can be expressed in exactly the same form as for 
the ID modulational case [|J, viz simply as a stabilizing 
damping according to 



where r(0,fc) « fcV3 — fc 2 /2, cf Appendix. This im- 
plies two things: the instability is suppressed by the 
incoherence for all wave numbers in the range [0, fc c ], 
where the cut-off wave number, fc c , is given by fc c = 
y3 — 4p| and secondly the range of instability decreases 
monotonously with increasing incoherence. However, 
this simple monotonously suppressing effect of the partial 
incoherence on the transverse modulational instability is 
not of a general nature. An indication of this was found 
in where numerical investigations were made using 
Gaussian as well as Lorentzian coherence spectra. Some- 
what counter-intuitively it was found that for the case of 
a Gaussian spectrum, increasing incoherence actually in- 
creased the range of modulationally unstable wave num- 
bers and increased the growth rate in part of the unstable 
region. Only for sufficiently strong incoherence did the 
unstable wavelength range start to shrink and the growth 
rate to decrease and to ultimately vanish. Thus, it seems 
that the properties of the transverse modulational insta- 
bility depend crucially on the form of the incoherence 
spectrum. That this indeed is so will be shown analyti- 
cally in the subsequent paragraph. 



IV. RESULTS FOR A GENERAL 
INCOHERENCE SPECTRUM 

In general, a complete analytical solution of Eqs.JBJ 
seems impossible to find. However, important informa- 
tion about the properties of the solution can be obtained 
by considering certain moments of the equations. For 
this purpose, we integrate the coupled equations for U 
and V over x and p x . This yields 



r((u)) + k Py ((v)) = o 

T((V))-kp y ((U))=G-((nR )), 



(17) 



where double brackets ((•)) denote integration over x and 
p x . The consistency condition for the real part of the 
perturbation can conveniently be expressed as 



{{U))d Pv = ((nS( Px ))). 



(18) 



Thus, solving for ((U)} from Eg. (|17|l and inserting this 
into Ea. (jl8() . we obtain the dispersion relation for the 
transverse instability of incoherent solitons in the form 



kp y G- 



r 2 + (k P y 



((nS( Px ))) 
((ni? » ' 



(19) 



T{p ,k) = r(0,fc)-fcp , 



(16) 



We underline that G_ is determined by the coherence 
properties of the soliton background solution, but that 
the parameter Q may depend on the coherence spec- 
trum. Nevertheless, the result expressed by Ea . (|19fl is 
completely general and is valid for arbitrary form of the 
coherence spectrum. We emphasize that Ea. (|19|) is of 
the same form as the dispersion relation for the modu- 
lational instability of apartially coherent, but homoge- 
neous, background, cf U 13 j i n which case the parameter 
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Q is easily determined to be Q = 1. On the other hand, 
for the transverse instability of a partially incoherent soli- 
ton stripe, the proper value of Q can not be easily found, 
although we may state that Q < 1. For the special case 
of a Lorentzian spectrum studied above, we can take one 
step further in Ea. (|19J) to obtain a dispersion relation 



(r + k Po ) 2 = (2Q - k 2 /2)k 2 /2 



(20) 



where, however, the Q factor still remains to be deter- 
mined. The analysis of the previous section and the result 
of the Appendix indicate that, Q ~ 3/4, independently 
of the degree of incoherence i.e. independently of pg. 
In order to pursue this line of analysis for general forms of 
incoherence spectra, we will assume weak partial incoher- 
ence in the sense that the incoherence spectrum is very 
narrow i. e. po -C k. The integral of Ea. l|19|) may then 
be evaluated approximately for any (well behaved) inco- 
herence spectrum G{p y ). This implies that the function 
F(p y ) = kpy/iT 2 + k 2 p 2 ) multiplying G_ in the integral 
can be expanded around the shifted wave numbers ±fc/2 
to yield 



F(Py)[G(Py + k/2) - G{ Py ~ k/2)]dp y 



k 2 fc 4 / fc 4 



2 

rmsi 



(21) 



where $ = T 2 + fc 4 /4 and we have defined p 2 ms = (x 2 ) as 
the rms-width of the spectrum 

/OC /POO 
x 2 G{x)dx I I G{x)dx. (22) 
-OC / J — OO 

The dispersion relation given by Ea. (|19fl then becomes 

r2 - fc2 («-T)-^( 3 -¥)^ ™ 

Since the incoherence is assumed weak, we will assume 
that the dispersion relation given by Eq. (|21(l may be sim- 
plified perturbatively by taking Q equal to its coherent 
value Q c and replacing <f> = T 2 + fc 4 /4 « k 2 Q c in the 
incoherently induced correction term. This yields 



k 2 

Q c - k - 



s fc2 i 



2 

rms 



(24) 



From this approximate expression for the growth rate, we 
can draw two important conclusions, valid for arbitrary 
(but narrow) incoherence spectra with finite rms-width: 
(i) the instability tends to be suppressed for all wave 
numbers in the range < k 2 < 3Q C , whereas in the re- 
gion 3Q C < k 2 < 4Qc, the growth rate is enhanced by the 
partial incoherence, (ii) the critical (non-zero) wave num- 
ber, fc c , at which the growth rate goes to zero, increases 
and is given by k 2 w 4(Q C + p 2 ms ) « 3 + 4p 2 ms . These 
analytical results agree well with what was obtained by 
Torres et al. Q using numerical computations. 



On the other hand, these results are in contradiction with 
the results obtained in the previous section for the case 
of a Lorentzian spectrum. There it was found that (i) 
the growth rate decreased for all wave numbers, (ii) the 
cut off wave number, k c , monotonously decreased with 
increasing incoherence. The explanation of this appar- 
ent contradiction is that the analysis of this section ex- 
cludes spectra, which, like the Lorentzian, do not have 
a finite rms-width. A direct implication of this result is 
that the effect of partial incoherence depends crucially 
on the form of the incoherence spectrum, even to the ex- 
tent that in some wavelength range the instability may 
even be enhanced by the incoherence. As demonstrated 
in |3j, for increasing incoherence, the range of unstable 
wave numbers first increases, but then eventually shrinks 
until finally the instability is completely quenched. This 
complete behavior is outside the range of validity of the 
perturbation analysis presented in the current section. 



V. CONCLUSION 

The present analysis has, in some detail, considered 
the effect of partial incoherence on the transverse mod- 
ulational instability of soliton stripes. We have shown 
that, for a Lorentzian form of the incoherence spectrum, 
the effect of partial incoherence on the transverse insta- 
bility agrees qualitatively with the corresponding result 
derived for the case of ID modulational instability; the 
growth rate decreases monotonously for increasing par- 
tial incoherence. However, the Lorentzian form is a very 
special case in the sense that although it has the nice 
property of being analytically integrable, it does not have 
a finite rms-width. Our analysis of general spectra with 
finite rms-widths shows quite a different qualitative be- 
havior of the growth rate for weak increasing incoherence. 
The growth rate is found to decrease for transverse wave 
numbers in the range < k < fc», but to increase in the 
complementary range fc* < k < fc c , where fc c is the cut 
off wave number of the instability and fc* is a characteris- 
tic transition wave number. In addition, it is found that 
fc c does in fact increase. These analytical results agree 
well with numerical simulations performed by as well 
as with previous analytical work of ours for the simpler 
case of the ID modulational instability, 



APPENDIX 

The dispersion relation for the transverse modulational 
instability cannot be determined analytically even in the 
coherent case and several different approximations have 
been presented, cf [3, [HA El]- We will here give a sim- 
ple, accurate and as far as we know new, approximation 
using a direct variational approach. Linearization of the 
2D coherent NLS equation, given in Eq.Q, around the 
stationary solution tp — sechccexp(iz/2) gives rise to two 
coupled equations for the real u(x), and imaginary v(x), 
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parts of the perturbed wave field. Inserting the assumed 
variations in y and z for the modulational perturbations 
(i. e. u, v oc exp(iky + Tz)), these equations become: 

Tu = L\v ; Tv = —L211 (A.l) 

where the operators L\ and L2 are self-adjoint and de- 
fined by 

f 1 d2 1/, , ,2n 1 2 

Ll =;2^ + 2 (1 + fc ) " SCCh *' (A.2) 

1/2 = 1/1 — 2 sech 2 x. 

can be reformulated as a variational problem 
corresponding to the Lagrangian Jz? = \vL\v — \uL2U — 
Tuv. An ansatz is made for the functions u and v as 
u = acf>2, v = /3</>i, where <px and <p2 are trial func- 
tions and a and /3 are the variational parameters. In- 
serting this ansatz into the variational integral, we find 



(J2f) = i« 2 (0i|Li|0i) - ajar^il^,) - i/? 2 (0 2 |i 2 |^ 2 ), 
where brackets (•) denote integration over x. The vari- 
ational equations with respect to a and (3 give rise to a 
linear system of equations for these parameters. A non- 
trivial solution of the system requires its determinant to 
vanish, giving the following dispersion relation 

2 _ (fo|£ 2 |fo)(fr|£i|0i) , A „, 
~ W^ 2 ' ( J 

With the intuitive choice of the trial functions as equal 
to the eigenfunctions of the operators L\ and L2, i. c. 
4>\ = sech a; and cf>2 = sech 2 a; respectively, the dispersion 
relation for the coherent case of the transvere instability 
becomes 

T 2 (0, k) « T 2 = fc 2 (3 - fc 2 )^ » ^(3 - k 2 ). (A.4) 
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